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The equations of motion for the molecular rotation are derived for vibrationally cold dimers that are polarized 
by off-resonant laser light. It is shown that, by eliminating electronic and vibrational degrees of freedom, a 
quantum master equation for the reduced rotational density operator can be obtained. The coherent rotational 
dynamics is caused by stimulated Raman transitions, whereas spontaneous Raman transitions lead to decoher- 
ence in the motion of the quantized angular momentum. As an example the molecular dynamics for the optical 
Kerr effect is chosen, revealing decoherence and heating of the molecular rotation. 
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I. INTRODUCTION 

Nonlinear optics, as commonly understood, describes the 
effects of nonlinear media on the optical field via suscepti- 
bilities, that are derived by a perturbative approach QJ. It is 
assumed in this context, that the medium can be treated as 
being stationary, which is well justified for thermalized sam- 
ples. The medium may be specifically composed of atoms, 
molecules or may be implemented as a solid-state device, such 
as a semiconductor The change of the quantum state of that 
physical system during the Ught-matter interaction is however 
largely neglected in the nonlinear optics literature, naturally 
since the focus lies on the properties of the light. Clearly the 
treated problem is a coupled one, where a description by sus- 
ceptibilities is the outcome of a decoupling of medium and 
light degrees of freedom that is only justified by certain re- 
strictive assumptions on the properties of the medium. A quite 
opposite viewpoint is to consider the dynamics of the atomic 
systems under the influence of the applied laser field 0]. In 
the latter case the scattered radiation is of minor importance, 
only its effect on the atomic dynamics by absorption and sub- 
sequent stimulated and spontaneous emissions of photons is 
taken into account. Thus the ancillary medium of nonlinear 
optics becomes the central issue. 

Focusing on the atomic or molecular dynamics however re- 
quires a treatment of physical conditions that are more gen- 
eral than those usually considered in nonlinear optics. For in- 
stance, if significant dynamics emerge, non-stationary initial 
quantum states and/or sufficiently strong coherent laser drives 
prevail. For atoms such situations have been reached in sev- 
eral implementations, the majority of them by providing low 
temperatures and atom-light coupling strengths comparable or 
larger than spontaneous decay rates. Prominent examples are 
atoms in high-Q microwave f3','4^,'5| or optical 01 cavities and 
single ions in radio-frequency traps 

Recent advances in trapping and cooling of atomic 0, [l3 
[HI d Q and molecular O El E O El gases, pro- 
gressively give access to the quantum-statistical properties 
and quantum dynamics of the manipulated gases. The ques- 
tion therefore arises as to how the quantum state of molecu- 



lar gases, for example, can be manipulated by nonlinear opti- 
cal interactions or how some of its properties can be revealed 
from the nonlinear optical spectra. For this purpose it is neces- 
sary to go beyond standard calculations of nonlinear optics, in 
order to treat not only the optical-field dynamics, but to focus 
on the dynamics of the molecular system. 

A first step towards such considerations is the study of 
the rotational dynamics of molecules that interact with off- 
resonant laser fields. Since dimers have no permanent dipole 
moment, they are infrared inactive and thus the rotational dy- 
namics is expected to be to large extent undamped. However, 
the polarizability of the molecule enables two-photon transi- 
tions and it will be shown here that these may in fact lead to 
decoherence in the rotational dynamics. 

Another aspect makes this issue particularly interesting. 
Consider ultracold molecules in a trap that is based on off- 
resonant optical fields, such as a dipole trap [15]. Assuming 
the molecules are cool enough to form a degenerate Bose gas, 
the lifetime of this peculiar quantum state is of major inter- 
est. As we will see here, off-resonant laser fields lead to a 
decoherence in the rotational quantum state, meaning that af- 
ter some time the initial rotational quantum state evolves into 
a statistical mixture of differing angular-momentum states. In 
that case one may expect that also the degenerate condensate 
state of molecules will be destroyed. 

In this paper we derive the basic equations of motion that 
describe the molecular rotational dynamics. We will focus on 
homonuclear, diatomic molecules interacting with both weak, 
off-resonant laser fields and electromagnetic vacuum modes. 
In such a configuration, besides the typical elastic Rayleigh 
scattering, stimulated and spontaneous Raman scattering pro- 
cesses occur that affect the molecular, ro-vibrational quantum 
state and may lead to decoherence. We derive a master equa- 
tion that includes these processes for studying rotational de- 
coherence, which may play an important role in the context of 
recent experiments with ultracold molecules flsll . The theory 
is applied to the optical Kerr effect, for which we solve the 
master equation by means of quantum trajectories. 

The paper is structured as follows: In SeclUthe physical 
situation under study is explained and the basic equations for 
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FIG. 1: Off-resonant laser excitation scheme of the molecule: \g) 
and |e) are electronic ground and excited states with respective inter- 
nuclear potentials, A is the detuning of the laser frequency ujq to the 
bare electronic transition frequency. 



the molecule-light interaction are given. Sec. Ulllis devoted 
to the adiabatic elimination of electronically excited states. In 
Sec. lIVI the trace over the vibrational degree of freedom is per- 
formed, leading to the reduced rotational master equation. It 
is solved in Sec.|V]for the case of an optical Kerr interaction 
and the damping and decoherence effects in the rotational de- 
gree of freedom are illustrated. Finally in Sec. lVIl a summary 
and conclusions are given. 



n. MOLECULE-LIGHT INTERACTION 

We consider a cold, homonuclear molecule interacting with 
a light pulse far detuned from possible resonances connected 
with electronic-state transitions of the type ^ ^E. More 
precisely, the molecule is assumed to be initially in its elec- 
tronic ground-state potential with a given distribution of pop- 
ulations among the ro- vibrational states. The laser field is then 
considered to be off-resonant, if all allowed ro-vibrational 
transitions from this distribution to excited electronic states 
are off -resonant, cf. Fig. [2 In this way two-photon Raman 
transitions dominate the dynamics and we consider in the fol- 
lowing stimulated as well as spontaneous Raman processes. 

Since electronic, vibrational and rotational transition fre- 
quencies are different from each other by orders of magnitude, 

IVeg :>UJ^-> B, (1) 

the dynamics reveals a separation of timescales. In Eq. Q ujeg 
is the bare transition frequency between the ground and near- 
est (dipole-allowed) excited electronic potential minima. The 
frequencies uj,^ and B are the typical vibrational frequency 



and the rotational constant, respectively. To avoid vibrational 
Raman transitions the spectral width Aw of the laser pulse 
obeys then the requirement 

Aw <c tj,., (2) 

and is typically of the order of rotational frequencies: Aw ~ 
B. Typically we may think of the laser field as being bichro- 
matic for addressing particular rotational transitions, though 
also more general spectra are possible. Note, that condi- 
tion (|2j is contrary to the situation in vibrational wavepacket 
creation by ultrashort laser pulses |2^ . 

We assume the molecule being located at a fixed posi- 
tion and omit its center-of-mass coordinate. Thus with re- 
spect to its center of mass a point-like molecule is considered. 
However, our derivations can in principle be extended to in- 
clude the quantum mechanical center-of-mass motion of the 
molecule. This extension would apply for describing molecu- 
lar beam deflecti on [[l9l I20I1 and would extend the purely co- 
herent treatment Bill by including spontaneous processes. 

We start with the general dynamical equations of a 
molecule interacting with a generic laser field and sponta- 
neously emitting photons into the vacuum electromagnetic 
reservoir The master equation for the density operator of elec- 
tronic, vibrational and rotational degrees of freedom q can be 
written as 

'Q = -l[H,g]+LQ. (3) 
n 

The time-dependent Hamiltonian reads 

H{i) = Hm + HuLit) , (4) 

where Hm is the free Hamiltonian of the molecule and 
Huhit) describes the interaction of the molecule with the ap- 
plied classical laser field. 

The Hamiltonian of the molecule can be written in the form 

Hm =^fi^g\9){g\ +^fii^e\e){e\ , (5) 
a e 

where \g) and |e) are short-hand notations for ro-vibrational 
energy eigenstates in the ground and excited electronic poten- 
tials, respectively, 

\9) = \9,Vg,3g,mg), \e) ^ \e,i^e,je,me) . (6) 

Our derivation applies to general molecular energy spec- 
tra. Only the ro-vibrational coupling needs to be neglected 
here. This coupling has been shown to lead to a decohering 
effect in the vibrational motion i22l I23I1 and thus can be ex- 
pected to produce decoherence also in the rotational degree 
of freedom. It appears, however, on a much larger timescale 
compared with the rotational frequencies and can be neglected 
if the molecule is and stays vibrationally cold during the in- 
teraction. Spontaneous Raman processes lead to a diffusion 
in the vibrational state, which may heat the vibrational de- 
gree of freedom. The diffusion rate, however, is small enough 
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for being neglected. In simplest approximation, for example, 
assuming a harmonic inter-nuclear vibration and taking into 
account only lowest-order rotational contributions, the eigen- 
frequencies could be given as 

= ^ly.gil^g + l) + Bgjg{jg + 1) , (7) 
UJe = UJeg+UJu.ei^e + ^)+Beje{je + i)- (8) 

Here io^^g/e are the vibrational frequencies in the ground and 
excited electronic potential, and i?g/e are the rotational con- 
stants in these electronic states. 

Turning back to Eq. @ the interaction of the molecule with 
the applied laser field E(t) reads 



where di is a vector component of the dipole moment and 



i7i„t(t) = -d.E(t) , 



(9) 



where d is the dipole operator of the single molecule. The 
laser field can be decomposed into slowly varying positive and 
negative frequency components E'^^)(t) by using the ansatz 



E(t) = E(+)(i) e 



-iujQt 



E(-)(t)e 



iujot 



(10) 



where ujq is the mid frequency of the laser pulse. Then the 
vector amplitudes E(^^(t) oscillate only at a maximum fre- 
quency given by the spectral width Auj, which is of the order 
of rotational frequencies. 

The incoherent part of the master equation Q describes 
spontaneous emissions of photons into the vacuum electro- 
magnetic reservoir and is given by the Liouville operator Cs- 
For the molecular energy levels being degenerate with respect 
to the rotational quantum number m the action of on the 
density operator reads ll24ll 



(15) 



For Eq. il5\ tiny modifications due to ro-vibrational frequen- 
cies have been neglected, since they are much smaller than the 
bare (optical) electronic transition frequency uieg- 

Note that the Liouville operator looks different from 
the usual cases considered in standard quantum-optical sys- 
tems, such as two-level atoms. This is only due to the ener- 
getic degeneracy with respect to the rotational quantum num- 
ber m. To properly account for this degeneracy we have to 
allow for subtle coherent effects: Starting from a ro-vibronic 
state |e, j/g, je, TOe), the spontaneous emission of a photon 
of frequency u! may leave the molecule in a coherent super- 
position of ro-vibronic quantum states \g, Vg^jg, rUg) with a 
unique value for Vg and jg but multiple values m £ [— Jg, jg] 
according to the selection rules. 

Given the dynamical equation (|3ji of the molecule interact- 
ing with the off -resonant laser pulse, in the following sections 
we will derive a master equation for the reduced rotational 
density operator ct, of the form 



n 



(16) 



The effective two-photon Hamiltonian Hn will contain the co- 
herent drive due to optically stimulated rotational Raman tran- 
sitions. It contains the polarizability tensor of the molecule 
and the classical laser field. The incoherent effects due to 
spontaneous Raman scattering will be described by the Li- 
ouvillean operator CsK{t)- 



-Cs ^ = ^ ^ ^ |a)(6| 5{uJab-i^a'b') 
i—1 ah a'b' 

X [{ama'){a'\e\h'){b'\^l\b) - \{a\^l^M'){a'\e\h'){b'\b) 



Ua\a'){a'\Q\b'){b'\lh\b) 



(11) 



The function 5{uj~uj') is a quasi delta function that ensures 
energy conservation www' for optical frequencies. 



5{lu) 



1 for |a;|<wo, 
else. 



(12) 



It selects resonant terms resulting from the performed secular 
approximation. Furthermore, the notation 



(13) 



denotes transition frequencies between the ro-vibronic molec- 
ular states \a) and \b) (a, b = g, e). 

The operators % reflect the spontaneous emission of a pho- 
ton due to molecular dipole moment pointing into direction i. 
That is, i labels the polarization of the spontaneously emitted 
photon. These operators are defined as 



7« 



K dj 



(i = l,2,3), 



(14) 



III. ELIMINATION OF EXCITED ELECTRONIC STATES 

For the interaction with an off-resonant laser pulse the 
molecule is hardly driven from the ground electronic poten- 
tial surface into the excited one. This applies, given that not 
too high vibrational levels are populated, which is the case 
for cold molecules as considered here. Thus the electronically 
excited states |e) = |e, I'e, je, JTie) can be adiabatically elim- 
inated to second-order in the molecule-light coupling. For 
this purpose we start with the coupled equations of motion of 
density-matrix elements, derived from the master equation Q 



Sag' 



' Qgg' 



-iujQtl 



^^li.geQee'llg'e' . (17) 

ee' 



? ee 



Qeg — 



'^'^egQeg H~ ^ ^ ^ ^eg' Qg' g ^ 
9' 



i ^ ^ Pee' ^e' f; ^ 2 ^ ^^9 



(18) 



(19) 
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and Qge = gig. Here gat are the matrix elements of the density 
operator in energy eigenstates 

gab = {a\g\b) (20) 
= {a,iyaja,ma\ g\b,iyb,jb,mb) {a,b^g,e), 

and 7i.ge are the matrix elements of the spontaneous-emission 
operators 7^: 

It.ge = {g\li\e) ■ (21) 

Furthermore, in Eqs jl8> and ( I19> the total decay rate F is 
defined as 

r^i:i:7:..7.,^-=v.5^^, (22) 

with d being the matrix element of the dipole moment of the 
considered electronic transition (cf. App. Finally, the 
time-dependent single-photon Rabi frequencies ^eg{t) are de- 
fined as 

ne,(t) = i(e|d|.g).EW(0 . (23) 



9gg' = 



The time dependent Rabi frequencies fl^g (t) oscillate much 
slower than the molecular vibrational frequencies since Aui <C 
iOjj gi^. They are slowly varying as compared to the oscillating 
exponentials in the integrals, that on average oscillate with the 
bare laser detuning defined by 

A = tJo - Weg . (27) 



They are slowly varying at frequencies given by the spectral 
width Aw of the laser pulse. 

We now formally integrate the equation of motion for the 
electronic coherences, Eq. ( I19> , from an initial time t' to a 
later time t, that obey the condition 

(24) 

and insert the obtained solutions into the equation of motion 
for the ground-state populations, Eq. \Y1\ . We consistently 
keep only terms up to second-order in the molecule-light inter- 
action. Defining the slowly varying density-matrix elements 

ea6W = ^ah(t)e-'""^*, (25) 

the approximation gab{t~T) « gab{t) can be performed for 
times r within the integration interval t — t' . From this pro- 
cedure the equation of motion for the electronic ground-state 
manifold is obtained. 



Thus fieg(i — r) w 0,e,g{t) for the time integrations and there- 
fore the Rabi frequencies can be taken outside the integral. 
After performing the remaining integrations and neglecting all 
remaining terms that oscillate with frequencies of the order of 
the detuning A or the vibrational frequency cu^ g/^, we obtain 



i ee' K. e" e' ^ 

•^0 



(c.c. and 5^5') 



(26) 



£'99' 



'Sag' 



VL*,{t)VL^g,{t) 
ggg"5u,v,, — - — (-lA-^r) + (c.c. and.g^.g') 



* nig{t)^,,g,{t) 



A2 



(28) 



Here Vg and Ve are the numbers of vibrational quanta in the 
ground and excited electronic potentials, respectively, and the 
Kronecker delta 5^^„ and 8^^^^, result as resonances from 
the rotating-wave approximation. This rotating-wave approx- 



imation with respect to vibrational frequencies holds here 
because we will consider from now on only density-matrix 
elements diagonal in the vibrational quantum numbers (in 
Eq. ( I28t Vg ~ Vgi). Those matrix elements do not freely os- 



5 



cillate with vibrational frequencies and since also the electric- 
field spectrum is narrower, they evolve on a slow timescale 
given by the Raman drive. Moreover, we have neglected small 
variations of the detuning due to ro-vibrational frequencies, 
since in the far off-resonant case considered here we may 
safely assume the condition 

A » Bg/, . (29) 

Considering a laser pulse that is well off-resonant to the 
electronic transition means also that the bare detuning is larger 
than the natural linewidth of the dipole transition by orders of 
magnitude, i.e., 

A > r . (30) 

Since for the case of a homonuclear molecule typically 
Bg/e ^ r, for the matter of consistency the original denom- 
inators in Eq. \2%\ have been further approximated by use of 

[A2 + (r/2)2]-iwA-2. 



Equation (I28> is still coupled to the electronically excited 
states Qee' ■ Thus, the manifold of density-matrix elements of 
electronically excited states is required. Performing the same 
procedure as for arriving at Eq. (I28> for the electronically ex- 
cited states, we obtain from Eqs ( fTsi and ( fT9t . 

Bee- = 2^ Qgg-5y,u^, ^5 ■ (31) 

gg' 

Again this result is accurate to second-order in the molecule- 
light interaction and is valid for matrix elements diagonal in 
the vibrational degree of freedom {v^ ~ v^i). This term rep- 
resents the excitation of the excited state from the electronic 
ground-states by absorption of a laser photon. 

Inserting OH into ( I28t and omitting terms of order higher 
than second in the molecule-light coupling, we arrive at 
the master equation for the consistently decoupled electronic 
ground-state manifold {Vg = Vgi) 



Qgg' = -i^gg'Qgg' + X! 

eg" 



Qgg"Sij^ij ,, 



A2 



(-iA-r/2) + (c.c. and ,9^,9') 



2^2^ 6'3"g"' Oi/g/zi/g/// ,2 l^,geli,g'e' 



(32) 



i ee' g"g"' 



This equation of motion describes the effects of spontaneous 
and stimulated Raman scattering processes on the dynamics of 
the ro-vibrational quantum state of the molecule. It represents 
a consistent approximation to second-order in the molecule- 
light coupling and preserves the trace of the density operator 
of the ground-state manifold: 

^Sgg^ XI i9,Vg,jg.'mg\g\g,Vg,jg,mg) =0 . (2,2,) 

Note, that the eigenfrequencies LOg, that appear in Eq. (I32> 
as transition frequencies ujggi = LOg—Ugi, only depend on the 
quantum numbers Vg and jg, i.e., 

UJg-^UJy^jg, (34) 

which is due to the energetic degeneracy with respect to the 
rotational quantum number vig. 

IV. ELIMINATION OF INTER-NUCLEAR VIBRATION 

Now the vibrational degree of freedom has to be eliminated 
to obtain the reduced dynamics of the molecular rotation. We 
intend to obtain an equation of motion for the reduced rota- 
tional density operator a, with matrix elements 

o-jm,j'm' = {jm\cF\j'm') , (35) 



by tracing the electronic ground-state density matrix Qggi over 
the vibrational degree of freedom, 

oc 

o-jm,/m' = X {9, fgj, m\g\g, Vg,j\ m') . (36) 

For notational convenience we have omitted here the 
electronic-ground-state subscript g for the rotational quantum 
numbers: jg, irig — > j, m. 

In Franck-Condon approximation the dipole operator is in- 
dependent of the inter-nuclear distance and therefore does not 
act in the Hilbert space of the vibrational degree of freedom. 
Then the Franck-Condon factors, being overlap integrals of 
vibrational wavefunctions in the ground and excited electronic 
potentials {vg\ve), will naturally emerge. Whereas, due to 
the narrow spectrum of the laser pulse, the stimulated Ra- 
man process only drives rotational transitions, spontaneous 
Raman scattering is accompanied by all possible vibrational 
Franck-Condon transitions. However, performing in Eq. (13 2> 
the trace over the vibrational quantum number Vg {vgi — Vg), it 
comes out that completeness relations of vibrational states can 
be used, so that combinations of Franck-Condon factors are 
summed over to result as unity. Thus it is possible to obtain 
an equation of motion for the rotation alone, that is no longer 
coupled to the vibrational dynamics. 

After having performed the sum over the number of vibra- 
tional quanta, the resulting master equation for the reduced 
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rotational density matrix can be written in the following Lind- 
blad form 



a ^ --[HM + HR{t),a] 



(37) 



J2 [m^sjit) ^islit)sM^^h<^slmit) 



The resulting time dependent Raman Hamiltonian HR {t) in 
Eq. ( I37K that describes the stimulated rotational Raman tran- 
sitions, reads 

Hnit) = E d |je,me)Oe,me|d.E(+)(t) . 

(38) 

It describes transitions from rotational states in the electronic 
ground state to those in the electronically excited ones and 
back to rotational states in the electronic ground state, con- 
nected with the absorption and stimulated emission of laser 
photons. 

Since the dipole operator acts due to its vector character 
also on the rotational degrees of freedom, we have to distin- 
guish between rotational states in the ground and excited elec- 
tronic states. Thus we have denoted this distinction in Eq. (I38> 
explicitely by the use of the subscript e. If no subscript is used 
the electronic ground state is assumed. 

The incoherent part of Eq. ( I37> describes the spontaneous 
Raman processes via the time dependent spontaneous Raman 
operator 

= rf.|je,me)(je,me|d.E(+)(i) , (39) 

with K being defined in Eq. il5i . From the derivation of 
Eq. (I37> it comes out that the Raman Hamiltonian ( I38> and 
the operator ( 13 9> are related by 



freedom. It reads 



(41) 



where the laser-coupling strength is determined by the two- 
photon Raman Rabi frequency. 



|d-E(+)p 



(42) 



(40) 



The Franck-Condon factors (z^li^e) for transitions between vi- 
brational state lu) and ji/g) of the ground and excited elec- 
tronic states, respectively, determine the strengths of spon- 
taneous Raman scatterings between vibrational states in the 
electronic ground-state potential. 

For vibrationally cold molecules, the internuclear poten- 
tials in the ground and excited electronic states may be ap- 
proximated by harmonic potentials, whose centers are shifted 
by Ax. Neglecting further the difference of vibrational fre- 
quencies in ground and excited electronic states, the Franck- 
Condon factors simplify to overlap integrals of shifted vibra- 
tional eigenstates. Starting from the vibrational ground state, 
the average vibrational excitation and variance can then be 
found from Eq. (gT} as i>{t)^gt and Av'^{t) = {gtf+gt{l + 
jjf), where p = iyy^r^Rr/A and rj^ /S.x^j2^ujy/h is the ra- 
tio of Ax to the spatial width of the vibrational ground-state, 
jjL being the reduced mass. 

From the condition uj^[v{t) + A.v{i)] <C A, which must 
be fulfilled in order that no resonant single-photon vibronic 
transitions occur, the valid dimensionless interaction times 
r = Bt, in units of rotational cycles B~^, can be obtained as 
r< [BA^/(rtj^l7R )]7?"2. In the weak-field regime, Or <i3, 
and for large detuning A, the first factor is rather large. Nu- 
merical values of 77 for alkali molecules are from 3.7 for ^Li2 
up to 8.6 for ^'^^Cs2. Thus, also the last factor 77"^ does not 
severly restrict the range of valid interaction times for our ap- 
proximations. 



This shows that the ratio of spontaneous and stimulated Ra- 
man transitions, inherent in the master equation ( I37> . is fixed 
by the ratio F/A and is independent of the laser inten- 
sity. Thus the laser field must be detuned by several natural 
linewidths from the electronic resonance to diminish sponta- 
neous processes compared to stimulated ones. 

A critical issue for the validity of our approach is given by 
the fact that we assume a far off -resonant laser field. Clearly 
spontaneous Raman scattering will lead to a diffusive broad- 
ening of the vibrational statistics that will slowly populate 
higher lying vibrational states. For our calculation to be valid, 
we do not have to assume that this process is negligible dur- 
ing the laser pulse duration. But we have to assume that this 
process will not populate higher lying vibrational states, that 
may eventually be subject to resonant transitions into an elec- 
tronically excited state. 

For estimating the increase of vibrational excitation during 
the laser interaction, a rate equation for the vibrational popu- 
lations can be derived by omitting the rotational degrees of 



V. ROTATIONAL DECOHERENCE IN OPTICAL KERR 
INTERACTIONS 

As an example for the decoherent effects of spontaneous 
Raman processes and as a first application of the master equa- 
tion derived here, let us consider a cw monochromatic laser 
field propagating along the z direction with linear polariza- 
tion in X direction. This configuration implements an optical 
Kerr interaction, where cr+ and (T_ polarized field components 
pick up different phase shifts depending on the rotational state 
of the molecule. The field is then subject to a rotation of its 
linear polarization while propagating through the molecular 
medium. Typically in the treatment of this optical effect the 
molecular medium is considered as being in a stationary (ther- 
mal) state, whereas here we focus on the actual dynamics of 
the molecule and do not consider the effects on the light field. 

Using the quantum- trajectory method, as described in 
App.IbI we have numerically calculated the dynamics of the 
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FIG. 2: Dynamics of the average angular- momentum component 
{Jy)- Initial quantum state is \j,9,(f>) with 6 = (p = it/2, and 
j = 2, i.e., the initial average angular momentum points into y di- 
rection. The laser interaction strength is SIr/B = 0.1 and detuning 
is A/T— 100. The dimensionless, scaled time is t = Bt; dotted line 
corresponds to F = 0. 



rotational degree of freedom. We start with an initial rota- 
tional quantum state that corresponds to an average angular 
momentum pointing in y direction. It is given by a coherent 
angular-momentum state, defined as 



E 

m=-j 



2j — m 
m 



X e 



cos(6i/2)]J'+"[sin(e'/2)]^-'" 
™*|j-m). (43) 



We take the quantization axis corresponding to the quantum 
number m as the z axis. For an average angular momentum 
pointing in y direction we then choose 9 = tt/2 and = 7r/2, 
which are the spherical angles of the average angular momen- 
tum. The quantum number j is chosen to be j — 2, as an 
example for a rotationally cold molecule that has been weakly 
excited. 

The numerically solved temporal evolution of the y compo- 
nent of the average angular momentum is shown in Fig.|2] The 
laser coupling strength is specified here by the two-photon Ra- 
man Rabi frequency fipj, cf. Eq. (I42> . which is chosen in units 
of the rotational constant B as ftn/B = 0.1. Moreover, the 
ratio of detuning to natural linewidth is chosen as A/F = 100, 
so that electronic transitions for a vibrationally cold molecule 
are widely suppressed. Nevertheless, rarely occurring sponta- 
neous Raman transitions are important for understanding the 
origin of decoherence effects, as demonstrated for a Raman- 
driven trapped ion ll26ll . 

Neglecting the spontaneous processes (F = 0) one obtains 
an oscillation of the average angular momentum as depicted 
by the dotted line in Fig. |2] In that case the average angular 
momentum has no projections on the x and z axes, so that a 
coherent one-dimensional oscillation is observed. The angular 
momentum changes from positive to negative y projections 
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FIG. 3: Time evolution of the fluctuations of the components of J for 
the same parameters as in Fig. 2(t — Bt; dotted line corresponds to 
r = 0). 



and vice versa via the z direction. This can be seen in its 
fluctuations, as shown for the dotted lines in Fig.|3] Whereas, 
the variance in x direction stays at a constant minimal value, 
a synchronized exchange of noise is observed in the variances 
of the angular-momentum projections in y and z directions. 

Taking into account the spontaneous Raman scattering of 
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FIG. 4: Heating of the rotational energy by spontaneous Raman pro- 
cesses in dependence on the scaled time r — Bt; parameters are the 
same as in Fig. 2. 



FIG. 5: Purity of the rotational state Tr(iT^) over scaled time r — Bt 
for the same parameters as in Fig. 2. 



the laser light that is detuned by 100 natural linewidth, co- 
herent processes are progressively suppressed, cf. solid lines 
in Figs 121 and |3] The oscillation of the average angular mo- 
mentum in y direction is damped, see Fig. |2] But it should 
be noted that no average components in x and z directions 
emerge throughout the evolution. Thus the overall coherent 
dynamics of the average angular momentum remains a one- 
dimensional oscillation, but it is eventually damped to zero. 

On the other hand, fluctuations in overall increase with pro- 
gressing time and therefore the average free rotational energy, 
being proportional to (J^), increases. This heating is caused 
by the spontaneous emissions of photons with random polar- 
izations and at random times. It is given by the sum of all 
three variances shown in Fig.|3] Whereas the oscillation depth 
of the noise in y and z directions gradually decreases with a 
compensating overall decrease of noise in y and increase of 
noise in z directions, there appears a monotonous increase of 
noise in the x direction. The latter is due to spontaneous Ra- 
man transitions that excite angular momentum in x direction, 
but with vanishing mean. Thus the overall rotational energy is 
subject to heating and increases as shown in Fig.|3] 

The decoherence effect of the spontaneous processes can 
be observed in the purity of the temporally evolving rotational 
density operator, see Fig.|5l Since any single (experimental) 
realization of the time evolution is interrupted by spontaneous 
emissions of photons with random polarizations and at ran- 
dom times, the ensemble average, describing the properties 
that can be expected from a typical experiment, gradually de- 
velops from an initially pure state into a statistical mixture. 



VI. SUMMARY AND CONCLUSIONS 

In summary we have derived a reduced master equation for 
the rotation of a small homonuclear molecule for the case of 



off -resonant laser excitation with pulse lengths larger than the 
period of the internuclear vibration. Several approximations 
have been justified and applied for obtaining the sought equa- 
tions of motion. First the electronically excited states have 
been adiabatically eliminated by assuming off-resonant laser 
fields for vibrationally cold molecules. Second the vibrational 
degree of freedom has been traced over, based on the Franck- 
Condon approximation. 

The obtained rotational master equation describes in its 
coherent part the interaction of the laser field with the 
orientation-dependent polarizability tensor of the molecule. 
This interaction is generated by stimulated Raman processes. 
The incoherent part, on the other hand, describes the sponta- 
neous Raman processes. The latter lead to a damping of the 
motion of the angular momentum of the molecule and, more 
importantly, cause substantial decoherence in the rotational 
quantum state. 

Our resulting master equation is of Lindblad form and can 
thus be efficiently solved by use of quantum trajectories. For 
the spontaneous rotational Raman transitions three different 
quantum jumps can occur, depending on the polarization of 
the spontaneously emitted photon. The importance of these 
spontaneous processes has been illustrated for the example of 
an optical Kerr interaction. Both damping of the rotational 
motion and quantum decoherence in the rotational state could 
be shown. 

Our results are relevant for ultracold molecules that are to 
be manipulated by off -resonant laser fields. An example is the 
interaction of the molecular rotational degree of freedom with 
the off -resonant laser fields in a dipole trap lll4lfl5ll . One may 
expect that even for trapping purposes decoherence appears in 
the rotational quantum state. Considering trapped molecules, 
our results hold for non-degenerate, non-interacting molecular 
gases. One may however guess, that the decoherence effects 
derived here will also emerge for degenerate molecular gases, 
where a full quantum-field theory applies. In that case one 
may expect a finite lifetime of, for example, molecular Bose- 



9 



Einstein condensates due to off-resonant laser interaction with which is independent of the ro-vibrational state, \iye,je, m-e), 
the molecular rotation, being at the lowest energetic scale. in the excited electronic potential. 
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APPENDIX A: TOTAL DECAY RATE 

For obtaining the total decay rate F, we start with the cal- 
culation of 



= EE(^l^ll-9)(5l7.|e')- (Al) 



APPENDIX B: QUANTUM TRAJECTORIES 



From Eqs ( I38> and ( I40> an effective non-Hermitean Hamil- 
tonian can be defined as 



iJcff (t) = + (^1 + i^j Hnit) . (Bl) 

Using this definition the reduced rotational master equa- 
tion (I37> can be rewritten as 



Hcsit) a ~ a Hlf,{t)\ + ^ ■ (B2) 



Using the definition (I14> we obtain 

ree'=«'EE(^l^ll5>(5Mde') 



(A2) 



t a 



For a (homonuclear) molecule having no permanent dipole 
moment, {g\di\g) = {e\di\e) = and thus the replacement, 

Y.\3)^9\^Y.\f^){9\ + Y.\^){e\^ i, (A3) 

9 g e 

can be performed, / being the identity operator, to obtain 

Fee' - K^e\d-d\e'). (A4) 

In Franck-Condon approximation the dipole operator does 
not depend on the internuclear coordinate. Moreover, d-d is 
a scalar that does not affect the rotational degrees of freedom. 
Therefore, we obtain 

Tee' = K, d Siy^iy'^ ^jeJe ^ ru ^ .m'^ ^ Tt^e.g', (A5) 

where d is the matrix element of the dipole moment of the 
considered electronic transition \g) ^ |e). The total decay 
rate F reads 



r = K^d^ 



(A6) 



This equation shows again the interpretation of Si{t) as being 
quantum-jump operators that are responsible for spontaneous 
Raman processes where a photon is emitted with polarization 
i. 

Therefore we may simulate quantum trajectories as realiza- 
tions of single experimental runs by sequential application of 
the non-unitary evolution operator 



U^s{t,t') = T exp 



(B3) 



where T denotes time ordering, and the jump operators Si{t). 
That is, starting from the normalized state |^(io)) at time to a 
quantum trajectory is given by 

\i){t\{{t, z)„, . . . , (t, z)i}) = ;7off(t, t„) (t„) . . . (B4) 

. . . S'», {t2) Ucs{t2,ti) S^, (tl) J/cff (tl, to) IV'(to)), 

where {t,i)k denote the time tk of a spontaneous Raman tran- 
sition and the polarization ik of the spontaneously emitted 
photon. 

Given the trajectory |V'(tn|{(t, i)„, . . . , (t, right after 

a spontaneous Raman process at time t„, the probability for a 
further spontaneous process to happen at a later time t is given 
by 



Pit) 



(V>(t|{(t, z)„, . . . , (t, z)i})|V>(t|{(t, . . . , (t, l)l})) 

{iitn\{{t,i)n, ■■■At, l)l})m„\{it, (t,l)l})) 

(?/'(t„|{(t,i)„,...,(t,i)i})|exp -^J^dTHnir) \ip{tn\{{t,i),, 



(B5) 



(V'(t„|{(t, i)„, . . . , (t, i)i})\^{t„\{{t, i)„, . . . , (t, i)i})) 
Drawing a random number R G (0, 1] the condition P{t) — R determines the time t = t„+i when the next spontaneous process 
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occurs. Then to decide which type i of spontaneous Raman process shall occur, the following probabilities are evaluated 

^ mtn+l\{it. l)n, ...At, l)l})\SUtn+l)S,{tn+l)mtn+l\{it, Ql})) ^ 1 2 3) (B6) 



I 

where According to the chosen polarization in+i the resulting tra- 

jectory is then 

mtn+l\{{t,t)n,---,{t,i)l})) (B7) 
= Ucs{tn+l,tn) \lp{t„\{{t,i)n, • ■ • , {t,i)i})), 

and W^ + VKj^ + VKj = 1. Drawing a second random number |V'(t„+i|{(t, i)„+i, . . . , (i, (B9) 

R' e (0, 1] the type of polarization in+i of the (71 + l)th spon- 
taneously emitted photon is chosen according to 

X ifO<R'<W^, 
i={ y ifW,<R'<W, + Wy, (B8) 

Z ifW^ + Wy<R'<l. 



= 5'j„+i(t„+i) \^Jj{tn+i\{{t,i)„, . . . , {t,i)i})). 
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